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An energy threshold of (220±10) eV was achieved at an efficiency of 50% with a four-channel 
ultra-low-energy germanium detector each with an active mass of 5 g. This provides a unique probe 
to WIMP dark matter with mass below 10 GeV. With a data acquisition live time of 0.338 kg-day 
at the Kuo-Sheng Laboratory, constraints on WIMPs in the galactic halo were derived. The limits 
improve over previous results on both spin-independent WIMP-nucleon and spin-dependent WIMP- 
neutron cross-sections for WIMP mass between 3—6 GeV. Sensitivities for full-scale experiments are 
projected. This detector technique makes the unexplored sub-keV energy window accessible for new 
neutrino and dark matter experiments. 

PACS numbers: 95.35.+d, 29.40.-n, 98.70.Vc 



There is compelling evidence from cosmological and 
astrophysical observations that about one quarter of the 
energy density of the universe can be attributed to Cold 
Dark Mattcr(CDM), whose nature and properties are 
still unknown[l(. Weakly Interacting Massive Particles 
(WIMP, denoted by x) are the leading candidates for 
CDM. There are intense experimental efforts^ to look 
for WIMPs through direct detection of nuclear recoils in 
%N^xN elastic scattering or in the studies of the possible 
products through \X annihilations. 

Supersymmetric(SUSY) particlesjl] are the leading 
WIMP candidates. The popular SUSY models prefer 
WIMP mass(m x ) in the range of ~100 GeV, though light 
neutralinos remain a possibility [4J. Most experimental 
programs optimize their design in the high-mass region 
and exhibit diminishing sensitivities for m x < 10 GeV, 
where an allowed region due to the annual modulation 
data of the DAMA experiment 0,0] — further reinforced 
by the first DAMA/LIBRA results Q — remains un- 
probed. Simple extensions of the Standard Model with 
a singlet scalar favors light WIMPs[8]. Detectors with 
sub-keV threshold are needed for probing this low-mass 
region and studying WIMPs bound in the solar system [§], 
and non-pointlike SUSY candidates like Q-balls[10]. This 
presents a formidable challenge to both detector technol- 
ogy and background control. Only the CRESST-I ex- 
periment has set limits llfl with sapphire(Al203)-based 



cryogenic detector at a threshold of 600 eV. 

The Kuo-Sheng(KS) Laboratory flij is located at 28 m 
from a 2.9 GW reactor core with an overburden of about 
30 meter-water-equivalence. Limits on neutrino mag- 
netic moments (/ij/)[T3| with a 1.06-kg germanium detec- 



tor(HPGe) at a threshold of 5 keV were reported[14|. 
These data also allowed the studies of reactor electron 
neutrinos [3| and reactor axionsflrjj]. A background level 
of ~ 1 event kg _1 keV _1 day _1 (cpkkd) at 20 keV, compa- 
rable with those of underground CDM experiments, was 
achieved. The current goal is to develop detectors with 
kg-scale target mass, 100 eV-range threshold and low- 
background specifications for the studies of WIMPs, \x v 
and neutrino- nucleus coherent scatterings [171]. 

Ultra-low-energy germanium detectors(ULEGe) is a 
matured technique for sub-keV soft X-rays measure- 
ments. They typically have modular mass of 5—10 g 
while detector arrays of up to 30 elements have been con- 
structed. Compared with AI2O3, Ge provides enhance- 
ment in xN spin- independent couplings (c^n) ^ ue to the 
A 2 dependence (H. fl"^ . where A is the mass number of the 
target isotopes. The isotope 73 Ge (natural isotopic abun- 
dance of 7.73%) comprises an unpaired neutron such that 
it can provide additional probe to the spin-dependent 
couplings of WIMPs with the neutrons(tr^). The nu- 
clear recoils from %N interactions in ULEGe only give 
rise to ~20% of the observable ionizations compared with 
electron recoils at the same energy. The suppression ra- 
tio is called the quenching factor(QF)fl9l|. For clarity, all 
ULEGe measurements discussed hereafter in this article 
are electron-equivalent-energy, unless otherwise stated. 

The ULEGe array consists of 4-element each having 
an active mass of 5 g 20j. Standard ultra- low-background 
specifications were adopted in its construction and choice 
of materials. It has identical external dimensions as 
the 1-kg HPGe of Ref. 14]. Apart from swapping be- 
tween the two detectors, data taking was performed with 
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FIG. 1: Schematic layout of the ULEGe with its anti- 
Compton detectors, inner shieldings and radon purge system. 
A 50-ton shielding structure encloses the detectors. 



all other hardware components, shieldings configuration, 
electronics and data acquisition (DAQ) systems [2l| kept 
identical. The schematic diagram of the experimental 
setup inside the shieldings is depicted in Figure [TJ 

The ULEGe signals were provided by built-in pulsed 
optical feedback pre-amplifiers, and were distributed to 
two spectroscopy amplifiers at 6 /zs(SAe) and 12 /is(SAi2) 
shaping times and with different amplification factors. 
Discriminator output of SAq defined the trigger condi- 
tions for DAQ. The threshold was set to about 4.3±0.2 
times the RMS fluctuations of the SAq signals above the 
pedestal. The DAQ rates for the ULEGe were about 
5 Hz, due mostly to electronic noise and agreed well with 
expectations [13, [23|. The SA 6 , SAi 2 , Anti-Compton 
Veto (ACV) and Cosmic-Ray Veto (CRV) signals were 
read out by 20 MHz Flash Analog-to-Digital Convertors. 
Random trigger(RT) events generated at 0.1 Hz and un- 
corrected with the rest of the system, as well as various 
system control parameters, were also recorded. 

Energy measurements were given by SA^ defined 
in the next paragraph. Figure [T^. shows an energy 
calibration spectrum due to external 55 Fe sources(5.90 
and 6.49 keV) together with X-rays from Ti(4.51 and 
4.93 keV), Ca(3.69 keV), and S(2.31 keV). Photons with 
energy lower than 2 keV were completely absorbed by the 
detector window. The RT-events provided the calibra- 
tion point at zero-energy. The RMS resolutions for the 
RT-events and 55 Fe peaks were about 55 eV and 78 eV, 
respectively. The calibration procedures were performed 
before and after the DAQ periods. Linearity was checked 
up to 60 keV with various 7-sources, and between zero 
and 2 keV with a precision pulse generator. The energy 
scale was accurate to < 20 eV, while deviations from lin- 
earity were < 1%. The electronic gain drifts, also moni- 
tored in situ by the pulse generator, were < 5%. A detec- 
tor hardware "noise-edge" of about 300 eV was achieved. 

Pulse shape discrimination (PSD) software was devised 
to differentiate physics events from those due to elec- 
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FIG. 2: (a) Top: Measured energy spectrum of the ULEGe 
with 55 Fe source together with X-ray peaks from various ma- 
terials. The black histogram represents events selected by 
PSD cuts. Deviations from the expected flat spectra con- 
tribute to PSD efficiencies, (b) Bottom: Scattered plots of 
the SAg versus SAj 2 signals, for both calibration and physics 
events. The PSD selection is shown. 



tronic noise, exploiting the correlations in both the en- 
ergy and timing information of the SAg and SA12 signals. 
Displayed in Figure [TJd is a scattered plot of the SAg and 
SAj 2 signals with the PSD cut superimposed, where the 
superscripts P/T denote partial/total integration of the 
pulses within (15,25) fxs and (-20,52) /is, respectively, rel- 
ative to the trigger instant (t=0). The noise events were 
suppressed. Calibration events and those from physics 
background were overlaid, indicating uniform response. 
Events selected by PSD but with CRV or ACV tags 
were subsequently rejected. The surviving events were 
ULEGe signals uncorrelated with other detector systems 
and could be WIMP candidates. 

Data were taken with the ULEGe at KS with differ- 
ent hardware and software configurations. They pro- 
vided important input on the background understanding 
and performance optimizations for future full-scale ex- 
periments. The data set with the best background and 
threshold has a DAQ live time of 0.338 kg-day. The anal- 
ysis results and the systematic effects at the two energy 
intervals which defined the sensitivities for m x below and 
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TABLE I: Summary of analysis results and their statistical and systematic errors in the two energy intervals which define the 
sensitivities at low and high WIMP masses. 



Energy Bin 198-241 eV 1.39-1.87 keV 



Raw Background Counts 


105212 


75 
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Trigger Efficiency (%) 




100 


DAQ Dead Time (%) 


11.0 ±0.1 


PSD — Cumulative Background Survival Fraction (%) 


0.008 


97 


Signal Efficiency (%) 


66 ±6 


100 


ACV — Cumulative Background Survival Fraction (%) 


0.0 


2.7 


Signal Efficiency (%) 


9£ 


.3 ±0.1 


CRV — Cumulative Background Survival Fraction (%) 


0.0 


0.0 


Signal Efficiency (%) 


91.4 ±0.1 


After-Cut Background Counts 








After-Cut Normalized Background Pates (kg - ^keV - ^day — ) 


± 27 2 (stat) ±|? (sys) 


± 1 g (stat) ± 0(sys) 


Quenching Factor 


0.200 ± 0.006 


0.244 ± 0.007 


Sampling in m x (GeV) 
a% (10- 39 cm 2 ) : 


5 


50 






Mean & Errors due to Background & QF Uncertainties 


0± ' 64 (Bkg)±0.01(QF) 


± 015 ,? (Bkg) ± 0.003(QF) 


Limit at 90% Confidence Level 


< 0.81 


< 0.20 


<rg° (10- 34 cm 2 ) : 






Mean & Errors due to Background & QF Uncertainties 


± 1£ >S (Bkg) ± 0.03(QF) 


0± ' 47 (Bkg) ± 0.01 (QF) 


Limit at 90% Confidence Level 


< 2.40 


< 0.59 
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FIG. 3: Selection efficiencies of the PSD cut, as derived from 
the 55 Fe-calibration and in situ data with ACV tags. Also 
shown are the best-fit la region and the trigger efficiency for 
physics events recorded by the DAQ system. 



above ~10 GeV are summarized in Table HJ 

The ULEGe data were taken in conjunction with a 
CsI(Tl) scintillator array[H| for the studies of neutrino- 
electron scattering. The combined DAQ rate was about 
30 Hz. The DAQ dead time and the CRV±ACV selec- 
tion efficiencies listed in TableUwere accurately measured 
using RT-events F4j . The maximum amplitude distribu- 
tions of the RT pedestals and physics events above the 
noise-edge of 300 eV were measured. The corresponding 
distributions for events <300 eV were evaluated by inter- 
polations to avoid biased sampling. The trigger efficien- 
cies depicted in Figure [3] correspond to the fractions of 
the distributions above the discriminator threshold level. 
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FIG. 4: The measured spectrum of ULEGe with 0.338 kg-day 
of data, after CRV, ACV and PSD selections. Background 
spectra of the CRESST-I experiment [10] and the HPGepJ] 
are overlaid for comparison. The expected nuclear recoil spec- 
tra for two cases of (m x , <t x n) are superimposed onto the spec- 
trum shown in linear scales in the inset. 



Events in coincidence with ACV-tags are mostly 
physics-induced. The fraction of these events surviving 
the PSD cuts was taken to be the PSD efficiency. This as- 
signment is conservative since the actual efficiency corre- 
sponds to the survival fraction of samples after electronic 
noise events in accidental coincidence were subtracted, 
and therefore should be higher. Alternatively, under the 
assumption that the 55 Fe calibration of Figure would 
give rise to physics events with a flat spectrum down to 
the lowest energy relevant to this analysis (<100 eV), 
the deviations of the PSD-selected events from a flat dis- 
tribution provided the second measurement. Consistent 
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FIG. 5: Exclusion plot of the spin-independent xN cross- 
section versus WIMP-mass, displaying the KS-ULEGe lim- 
its and those defining the current boundaries [HI, 28]. The 
DAMA-allowed regions^ are superimposed. The striped re- 
gion is that favored by SUSY models [J]. Projected sensitivi- 
ties of full-scale experiments are indicated as dotted lines. 



results were obtained with both approaches, as depicted 
in Figure [3l The larger uncertainties of the first method 
are due to the limited statistics from only the in situ ACV 
samples. The efficiencies and their uncertainties adopted 
for analysis were derived from a best-fit on the combined 
data set. A threshold of (220±10) eV was achieved with 
a PSD efficiency of 50%. 

The ULEGe spectrum normalized in cpkkd unit af- 
ter the CRV, ACV and PSD selections is displayed in 
Figure [H showing comparable background as CRESST- 
l[llj. Listed in Table [J are the normalized background 
rates, indicating that statistical uncertainties dominate 
over the systematic effects. The formalisms followed 
those of Ref. [l8| using standard nuclear form factors, 
a galactic rotational velocity of 230 kms -1 and a local 
WIMP density of 0.3 GeV cm -3 with Maxwellian veloc- 
ity distribution. No subtraction of background profiles 
was made such that the WIMP signals cannot be larger 
than the observed event rates. The unbinned optimal 
interval method as formulated in Ref. [25| and widely 
used by current CDM experiments was adopted to de- 
rive the upper limits for the possible %N event rates. 
By comparing the observed background in different en- 
ergy intervals with the expected number of events due to 
%N recoils for each m x , the optimal intervals producing 
the most stringent limits to crSj and cr™ 1 were selected. 
Corrections due to QF, detector resolution and various 



efficiency factors were incorporated. The energy depen- 
dence of QF in Ge was evaluated with the TRIM software 
package |26|. The uncertainties were taken to be their dif- 
ference with the statistical best-fit values of the available 
data[27| from 254 eV to 200 keV nuclear recoil energy. 



Exclusion plots on both (m x ,<r^) and 



D ) 



planes at 90% confidence level for galactically-bound 
WIMPs were then derived, as depicted in Figures [5] and 
[5k., respectively. The DAMA-allowed regions 0, Q and 



the current exclusion boundaries [111, |28j are displayed. 
The model-independent approach of Refs. [29|, [3(j were 
adopted to extract limits on the spin-dependent cross- 
sections. Consistent results were obtained when different 
73 Ge nuclear physics matrix elements [3l| were adopted as 
input. The constraints on the effective axial four-fermion 
X~ proton and x-neutron spin-dependent couplings [30l| at 
m x =5 GeV are displayed in Figure [J5Jd. The param- 
eter space probed by the 73 Ge in ULEGe is comple- 
mentary to that of the CRESST-I experiment where 
the 27 A1 target is made up of an unpaired proton in- 
stead. New limits were set by the KS-ULEGe data 
in both cr^jj and ct x ° for m x ~3— 6 GeV. The remain- 
ing DAMA low-m x allowed regions in both interactions 
were probed and excluded. The observable nuclear re- 



0.5 x 10" 39 cm 2 (allowed) 



coils at m x =5 GeV and 0" x n = 
and 1.5 x 10~ 39 cm 2 (excluded) are superimposed with 
the measured spectrum in the inset of Figure 2] for il- 
lustrations. It is expected that recent data from the 
COUPPjUj] experiment can place further constraints in 
the spin-dependent plots of Figures [6^&b. 

This work extends the bounds on WIMPs by mak- 
ing measurements in a new observable window of 
100 eV— 1 keV in a low-background environment. Under- 
standing and suppression of background at this sub-keV 
region is crucial for further improvement. Measurements 
are conducted with the ULEGe at an underground lab- 
oratory. There are recent advances in "Point-Contact" 



Ge detector [33[ which offer potentials of scaling- up the 
detector mass to the kg-range. Preliminary results in 
dark matter searches were recently reported [34j. The 
mass-normalized external background will be reduced in 
massive detectors due to self-attenuation [17]. Further 
reduction in threshold may be possible with improved 
junction field-effect-transistors and by correlating signals 
from both electrodes. The potential reach of full-scale 
experiments with 1 kg-year of data and a benchmark 
background level of 1 cpkkd is illustrated in Figures [5] 
and EH- Such experimental programs are complementary 
to the many current efforts on CDM direct searches. 

We are grateful to the KIMS Collaboration and the 
authors of Ref. [22j for inspiring comments. This work 
is supported by the Academia Sinica Pilot Project 2004- 
06, Theme Project 2008-10, contracts 95-2119-M-001-028 
and 96-2119-M-001-005 from the National Science Coun- 
cil, Taiwan and 10620140100 from the National Natural 
Science Foundation, China. 
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FIG. 6: (a) Left: Exclusion plot of the spin-dependent x-neutron cross-section versus WIMP-mass. Same conventions as those 
in Figure [5] are used, (b) Right: Constraints at m x = 5 GeV on the effective axial four-fermion x-P r °ton (a p ) and x-neutron 
(tin) spin-dependent couplings, in units of 2\/2Gf following the formulation of Ref. [30l ]. The parameter space within the bands 
of the corresponding experiments are allowed. The shaded area at the origin is the combined allowed region. 
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Appendix I : 

Trigger Efficiency at Threshold — 
Reply to arXiv:0806.1341k r2 

H.B. Li, S.T. Lin, S.K. Lin, A.K. Soma, H.T. Wong* 

(TEXONO Collaboration) 
* Corresponding Author: htwong@phys.sinica.edu.tw 

The authors in a recent paper [l[ raised questions on 
our estimates of trigger efficiency in the evaluation of 
the constraints on WIMP Dark Matter in the low-mass 
(< fO GeV) domain with a Ultra-Low Energy Germa- 
nium (ULEGe) detector at a threshold of 220 eV 0. 
The discrepancy originates from some misleading termi- 
nology in Ref. [2| . We address the issue in details in this 
reply. We show how the trigger efficiencies are derived 
and demonstrate that the results of Ref. Q are correct. 

I. CONCEPT OF EXPERIMENT 

Details of the experimental setup and data analysis, 
as well as the definitions of notations, can be referred to 
Ref. [2| . There are two categories of events relevant to the 
present discussion: physics signals (PHY) and electronic 
noise (ELE). The PHY events are due to actual energy 
depositions at the ULEGe by gammas, neutrons, neutri- 
nos, WIMPs and other radiations, while ELE events are 
induced by the various stages in the readout electronics. 
The majority of the events above the noise edge of 300 eV 
are PHY-events. A particularly good method to extract 
a clean sample of PHY-events below the noise edge for 
further studies is to require that the ULEGe signals are 
in coincidence with the Anti-Compton Veto (ACV) tag. 
However, PHY-events due to neutrinos or dark matter 
interactions would have to be in anti-coincidence with 
ACV. A major goal of the experiment is to lower the 
threshold by suppressing ELE while keeping PHY-events 
in some substantially large and known fraction. 

In order for PHY-events to be included in the final 
spectra where physics is extracted, they have to survive 
three procedures, the efficiency of each of which must 
be known. These efficiency factors were summarized in 
Table 1 of Ref. 0: 

1. Trigger efficiency (eTrig) — that PHY would pro- 
duce a trigger signal to the data acquisition (DAQ) 
system. This efficiency depends on the energy of 
PHY. 

2. DAQ efficiency (6daq) — that the trigger signal 
would actually activate the DAQ system resulting 
in a complete event recorded in the computer. This 
efficiency is independent of the energy of PHY. 

3. Analysis efficiency (eAniy) — that offline software 
procedures would retain PHY and suppress ELE. 



Some PHY would be rejected in the process leading 
to an efficiency factor which is energy-dependent. 

Among these three efficiency factors, the derivations of 
£daq and eAniy were discussed in Ref. [2J . The cdaq was 
evaluated accurately by "Random Trigger" (RT) events 
to be 89% - fraction of the RT events actually recorded in 
the computer, (cited in Table 1 as "DAQ Dead Time" of 
1 1 .0%.) Same procedures were used in our earlier work 
on neutrino magnetic moments 0]. The evaluation of 
£Aniy by two different approaches was discussed in the 
text and the results were shown in Figure 3. Further 
elaborations are made in Ref. 0]. The rest of this reply 
would focus on eTrig- 



II. DISCREPANCIES AND ORIGINS 

Ref. *!!] challenged our values of £daq ~ 89%, stating 
that the "discriminator threshold of 20 eV" and "RMS 
resolution of pedestal (that is, RT events) of 55 eV" 
would imply a DAQ rate (R) of ~20k Hz, based on a 
known relation [5J] 



R ~h exp[ ~4^ ] 



(i) 



where r is the shaping time (=6 [is for the trigger signal 
SAe), <7 is the RMS of the pedestal noise fluctuations 
and d is the threshold level above the pedestal. The 
calculated rate is much higher than the actual DAQ rate 
of ~5 Hz for the ULEGe in the actual measurement. 

This discrepancy originates from the fact that it is in- 
correct to use these two energy values together in Eq. [T] 
In the experiment, the "energy" of an event is defined 
by integrating (the "Q-mode") the SA12 signal within 
time intervals (-20,51.2) /j,s where t=0 denotes the trig- 
ger instant. That is, the energy is measured through the 
summation of 71.2X20=1424 FADC numbers. Under this 
definition, the pedestal (RT) events have an RMS resolu- 
tion of 55 eV, after appropriate calibration. However, the 
triggering was done by comparing the SAg signal with a 
pre-set discriminator level. That is, the relevant quantity 
in this operation is the amplitude of the pulse (the "A- 
mode") — and in fact of a different pulse (SAg instead of 
SA12, but this is a minor point). The RMS-resolution of 
RT events in Q-mode is not related to the trigger config- 
uration in A-mode and therefore should not be taken as 
the a of Eq.\l\which describes the trigger rate. It is an ex- 
perimental optimization frequently used (also discussed 
in Ref. H) — that the energy definitions are based on 
averaging over long periods, while the triggering schemes 
make use of the instantaneous response of the signal. 
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FIG. 1: Typical SAe pulses due to (a) Top: a random trigger 
event and (b) Bottom: a physics event at 139 eV in coinci- 
dence with anti-Compton signal. The pedestal mean and the 
discriminator threshold were denoted by dotted and dashed 
lines, respectively. 



III. EVALUATION OF TRIGGER EFFICIENCY 

Having identified the source of discrepancy which led 
to the incorrect conclusions in Ref. we present the 
evaluation of fTrig in what follows. Displayed in Fig- 
ures [IJi&b are typical SAg signals for RT and PHY 
events, respectively. The amplitude is presented in 
FADC unit (2 V is equivalent to 256 FADC unit). The 
PHY-event was measured to be 139 eV in the Q-mode, 
and in coincidence with ACV. Superimposed is the dis- 
criminator level (38±2 mV in hardware unit, uncer- 
tainties mostly from instabilities over time) which is 
4.8 FADC unit above the mean value of the pedestal level. 
It can be seen that the RT event is below threshold while 
the PHY event is above threshold by ~ 1.5 FADC unit, 
thereby provided a trigger to the DAQ system with good 
margin. 

Histograms of large samples of events like those of Fig- 
ures QJ,&b is presented in Figure^ where the horizontal 
axis is in FADC unit. The "Random Noise" histogram 
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FIG. 2: The distributions of noise fluctuation of RT-events as 
well as of the maximum amplitudes of PHY-events in various 
energy bins. The discriminator threshold level is also shown. 



corresponds to the distribution of the amplitudes of every 
time-bin of the RT-cvcnts. The RMS is 1.1 FADC unit 
and represents the noise fluctuation of the SAg signal. 
The discriminator threshold level is also shown. Accord- 
ingly, one can equate r=6 /is, d=4.8±0.3 and <r=l.l in 
Eq. [U giving R ~ 3.1 ~_2% Hz. The measured rate aver- 
aged over the entire DAQ periods of about 17 days and 
for four channels together is ^5 Hz. This is in good 
agreement. 

The various distributions show the maximum ampli- 
tude of PHY-events at various energy ranges measured 
by the Q-mode. All the events above 100 eV (the rel- 
evant range for subsequent analysis) exhibit at least 
1 FADC Unit of margin above threshold. The trigger 
efficiencies were then derived using the maximum am- 
plitude distributions of the RT events at E=0 and the 
PHY-events between 300 eV to 1000 eV. The mean and 
RMS for the E=0— 300 eV regions were evaluated by in- 
terpolation (rather than from actual data) to avoid bi- 
ased sampling. The results are displayed in Figure [3] 
We note that the energy range that provides the most 
severe constraints to the dark matter analysis is that of 
200-250 eV, where the trigger efficiency is close to unity. 

Using PHY-events as well as data taken with precision 
test pulser, it can be shown that energy measurements 
with Q-mode (integration) and A-mode (maximum am- 
plitude) are both valid and equivalent for pulses which 
are large compared to the noise fluctuations. Moreover, 
the Q-mode measurements are well-behaved and linear 
all the way towards the pedestal zero-energy level. On 
the contrary, the A-mode measurements become inac- 
curate and the calibration is non-linear as the energy 
decreases and approaches zero. Although the threshold 
is well-defined in amplitude (4-8±0.3 FADC unit), the 
statement "discriminator threshold at 20 eV" in Ref. 0/ 
has been misleading and too simplistic. It is incorrect to 
equate the threshold to a single number in Q-mode as its 
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1 Threshold = 4.25 ± 0.22 Noise-RMS 




0.25 



Energy (keV) 

FIG. 3: Trigger efficiencies as function of energy as measured 
by the Q-mode. 



energy scale. That is, the trigger efficiency versus energy 
curve of Figure [3] is not a step-function. For instance, 
£Tri g = 50 ± 30% corresponds to an energy range of about 
80±50 eV. This is the root of the misunderstanding and 
is corrected in the revised text. 



IV. COMPARISON OF KS AND Y2L SPECTRA 

Another point noted in Ref. [l| is the comparison of 
spectra in Ref. 0] with those taken at the Y2L under- 
ground laboratory in South Korea .6] . Such comparisons 
would not be appropriate. Different detectors were used 
even though the specifications are similar. The electronic 
modules and noise sources (which are affected by many 
ambient conditions) are not identical. The most impor- 
tant difference is that we did not make an attempt to 
use pulse shape differentiation (PSD) techniques to sup- 
press ELE-events with Y2L data. The spectra shown 
with Y2L data are always without PSD cuts, such that 
the raw "background" — dominated by ELE-events — 
below the noise edge are extremely high. DAQ rates in 
both cases are similar, at the range 1-10 Hz. Displayed 
in Figure S] are the raw spectra for both KS and Y2L, as 
well as that after PSD suppression for KS. It can be seen 
that the raw spectra are comparable. The residual differ- 
ences are due to ambient electronic noise and radioactive 
background conditions. 

We are indebted to Profs. F. Avignone, P. Barbeau 
and J. Collar for careful reading and critical comments 
on Ref. |2|- Their input not only expose the parts of 
the original text where the experimental procedures were 
not adequately explained, but also stimulate us to more 
thorough thoughts on the issue. 




Measured Energy (keV) 



FIG. 4: The raw energy spectra of Y2L and KS, as well as 
that after PSD suppression for KS. 
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Appendix II: 

Selection Cuts and Efficiencies at 

Threshol d — 

Reply to larXiv:0806.1341k f3 

H.B. Li, S.T. Lin, S.K. Lin, A.K. Soma, H.T. Wong* 

(TEXONO Collaboration) 
* Corresponding Author: htwong@phys.sinica.edu.tw 

Additional questions were raised in a recent paper [l[ 
on the experimental concepts and systematic issues of our 
recent results on WIMP Dark Matter searches in the low- 
mass (< 10 GeV) domain with a Ultra-Low Energy Ger- 
manium (ULEGe) detector at a threshold of 220 eV 0. 
We provide clarifications and justifications on these is- 
sues, and conclude that there are no flaws in our proce- 
dures. 

I. COMMENTS AND REPLIES 

In a previous version(V2) Q, the authors raised ques- 
tions on our DAQ dead time, trigger efficiencies and ap- 
parent inconsistencies between the Kuo-Sheng Reactor 
Laboratory (KS) and the Yang- Yang Underground Lab- 
oratory (Y2L) data. In an earlier reply [H, these ques- 
tions were addressed in details. No further comments 
were made along these lines. 

A new version(V3) of the comments [H was subse- 
quently posted, where questions on the experimental con- 
cepts and various systematic issues were raised. We re- 
sponded on these comments in this article. Section II 
deals with the basic concepts and performance of the 
PSD cuts and efficiencies, while Section III elaborates 
on the reasons behind the various choices made in the 
experiment. Details of the experimental setup and data 
analysis, as well as the definitions of notations, can be 
referred to Ref . 0] • 

II. SELECTION CUTS AND EFFICIENCIES 
A. General Comments 

The objectives of applying selection criteria ("cuts") 
on experimental data is to reject undesirable "back- 
ground" events and increase the fraction of "signal" 
events in the data sample. The events surviving the cuts 
need not be all signals, only that usually, the signal-to- 
background ratios are enhanced by the selection. 

There are much freedom in the choice of these cuts. 
There may be effective or ineffective cuts, but in general, 
all cuts are valid. The variables ( "figure of merits" FoM) 
on which the cuts are applied are defined to optimize 
the performance. The FoMs can be mathematical con- 
structions and need not correspond to, or be linear to, 
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FIG. 1: Scattered plots of the SAg versus SAJ2 signals, for 
both calibration and physics events before ACV and CRV 
cuts. The PSD selection is shown. 

certain physical variables. Once the events are selected 
by the cuts, the physical parameters (like energy) can be 
derived for these events with different algorithms which 
are themselves constructed to give the best resolutions 
for these quantities. 

When a set of cuts are applied to experimental data, 
the corresponding "selection efficiencies" must be eval- 
uated. These are the probabilities that potential signal 
events that will survive these cuts. The goals of apply- 
ing cuts are to suppress background events as much as 
possible while keeping the selection efficiencies as large 
as possible. 



B. Specific: This Analysis 

In this particular analysis [2|, the ACV and CRV cuts 
are straightforward. They are signals from detector com- 
ponents other than the ULEGe target, and their efficien- 
cies were evaluated by random trigger (RT) events. Their 
performance are summarized in Table 1 of Ref. • 

It is the PSD-cut which was questioned by Ref. 
This cut was applied to the variables "SA 6 " and "SA 12 ", 
as displayed again in Figure [TJ Among them, SA12 also 
provides an energy measurement, while SAg is a math- 
ematical construction to optimize the performance of 
signal-vs-background differentiation. In this case, "back- 
ground" corresponds to electronic noise events below the 
hardware noise edge. 

Motivations for the choice of : The two pulses S 
and SA12 are correlated for physics signals. The correla- 
tions arc different and less strong for noise-events. The 
conceptual idea behind the PSD cut is that, given the 
energy of an event is known through measurements by 
SA12, positive fluctuations of the SA6-pulses in physics- 
induced events at a particular time interval and ampli- 
tude ranges can be expected. This is probed by the 
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PSD cut. The selected time interval is optimized for 
energies near threshold (< 300 eV), such that this in- 
terval no longer corresponds to the amplitude-peaks of 
SA 6 for events at higher energy. Therefore, SA 6 being 
non-linear with energy does not jeopardize the validity 
of the cut. The reasons behind the non-linear behaviour 
is well-understood — and in fact intentional. 

Is the PSD cut arbitrary?. NO. The application of the 
PSD cut is based on genuine understanding on the be- 
haviour of the detector hardware. Such PSD techniques 

— correlations of two signals at two shaping times and, 
alternatively, of the full and partial integrations of the 
signals — are well-established ones at higher energy in 
the case of a/7 separation in many detector systems. 

Displayed in Figure [2] are the survival fraction (/) of 
events at E=200-300 eV with an ACV (Anti-Compton) 
tag versus the relative timing between the ACV signals 
and the ULEGe triggers. Overlaid are the actual coin- 
cidence time interval between the ACV and the ULEGe 
systems determined independently from hardware tim- 
ing. The coincidence window is about 3 /is. This is 
defined by (a) the 6 (is shaping time output SAg from 
ULEGe which provided the trigger timing; and (b) the in- 
creased time-jitters at low energy. The ACV and ULEGe 
signals outside this range are accidentals and uncorre- 
cted. 

It can be seen that ONLY ULEGe events in correct 
coincidence with the ACV-tags give a substantial value 
of /. The survival fraction for events without ACV-tags 
(denoted by the data point at negative time in Figure [2|) 

— predominantly due to electronic noise at this energy — 
is (1.7±0.3) x 10~ 4 . This proves that the PSD cut is suc- 
cessfully devised and is performing its intended functions 
of "suppressing electronic noise events and selecting the 
physics-induced events" . 

Selection Efficiencies: The PSD-efficiency for physics- 
induced events (epsd) is related to the measured survival 
fraction / by 

_ (epsp * P + }n * N) 
} ~ P + N 

where P and N are, respectively, the numbers of physics- 
induced and noise events in correct coincidence with 
ACV, while Jn is the survival fraction for events trig- 
gered by electronic noise. In this experiment, / and (P + 
N) are measured quantities, Jn ~ 10~ 4 at 200—300 eV 
and even less at lower energy, while the average N can 
be evaluated from the non-coincidence samples where 
P = 0. It can be seen that, in general, £psd > / so 
long as there are finite fraction of noise events in the 
ACV-tagged sample. 

Statistics are limited in the present analysis since only 
in situ data were used. The number of noise events in the 
coincidence ACV-tagged sample is finite but has large un- 
certainties. Accordingly, the more appropriate approach 
is to take zero noise-background (N — 0), giving rise to 
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FIG. 2: The survival fraction of events at E=200-300 eV 
with an ACV (Anti-Compton) tag versus the relative timing 
between the ACV signals and the ULEGe triggers. Overlaid 
are the actual coincidence time interval between the ACV 
and the ULEGe systems derived independently from hardware 
timing. The data point at negative time is due to events 
without ACV-tags, and corresponds to (1.7 ±0.3) x 10~ 4 . 

the assignment of epsD = /■ The subsequent upper lim- 
its derived would therefore be conservative and less con- 
straining ones. The variations of epsD with energy under 
this assignment are displayed in Figure 3 of Ref . 0] • 

III. EXPERIMENT'S CHOICES 

Inevitably, an experiment has to choose various tools 
from the available pool in the course of its analysis. We 
explain why such choices were made for the cases raised 
in Ref. [1], and illustrate the sensitivities to the physics 
results if alternative schemes would had been selected 
instead. 

A. Fe spectrum for Selection Efficiency 

It was explicitly stated in Ref. [2J that the flat 55 Fe 
spectrum at low energy is an assumption. Under this 
assumption, the PSD selection efficiencies were derived. 
The results agree well with those obtained by the more 
rigorous approach with ACV-tagged events. The good 
agreement suggests validity of the assumption at the 
present level of accuracy. Accordingly, it is justified to 
add these results to further constrain (to reduce uncer- 
tainties of) our knowledge of epsD- This is the approach 
adopted in Ref. 0]. 

If the data from the 55 Fe spectrum would be ig- 
nored altogether and the PSD efficiencies would be de- 
rived exclusively from the ACV-tagged events, epsD 
at E=200-300 eV would move from (0.65±0.06) to 
(0.61±0.08). The 50% efficiency line would correspond to 
224 eV instead of 216 eV, while the arL limits (in units 
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FIG. 3: A compilation of all quenching factor (QF) mea- 
surements on germanium, with calculations from the TRIM 
software [6] as well as by the Lindhard model [7| under two 
parametrizations (k=0.20 and 0.15) overlaid. 

of 1CP 39 cm 2 throughout in this Section) at m x = 5 GeV 
would increase (become less constraining) from 0.81 to 
0.88. 

B. Quenching Factor 



m x = 5 GeV 



C. Constructing Exclusion Plots 

The unbinned "optimal interval method" as formu- 
lated in Ref. [H was adopted to derive the exclusion lim- 
its. The unbinned formalism allows the use of all avail- 
able information in the background spectra and was used 
in other CDM experiments like CDMS and XENON. NO 
background profile was assumed or subtracted, which is 
also a conservative approach. The sensitivities at low m x 
under this scheme are driven by the absence of counts be- 
tween 198 eV and 241 eV. 

An alternative method would be to place the back- 
ground events in different energy bins and follow the for- 
malism of Ref. 9]. For instance, choosing 50-eV bins 
for E>100 eV (thereby deliberately filling the hole at 
200-250 eV), the <r^ limit at m x = 5 GeV would in- 
crease (become less constraining) from 0.81 to 1.20. This 
reduction in sensitivities is expected since data binning 
involves loss of information. 

We conclude that our choices in these three aspects 
of the experiment are justified. The sensitivities of the 
physics results (exclusion upper limits) are dominated by 
the statistical uncertainties of the background spectra. 
The potential effects on them are minor if alternative 
schemes would have be chosen instead. 



A compilation of all quenching factor (QF) measure- 
ments on germanium is given in Figure [3] Overlaid are 
calculations from the TRIM software Q as well as by the 
Lindhard model Q under two parametrizations (k=0.20 
and 0.157). Both schemes have been adopted in various 
CDM experiments. It can be seen that the TRIM re- 
sults explain well the QF measurements at both low and 
high energy. Accordingly, we chose to use this scheme in 
our analysis. The QF values are less than those evaluated 
with the Lindhard (k=0.20) model, and hence would give 
rise to more conservative results. 

If Lindhard (k=0.20) would be used, the QF at 1 keV 
recoil energy will be increased from 0.20 to 0.21. The 
QF uncertainty estimations of 0.006 in Ref. @] can ac- 
count for this deviation. This alternative choice will 
only have minor effects on the exclusion limits, decreas- 
ing it (becoming more constraining) from 0.81 to 0.80 at 
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